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QCD in 1+1 dimensions is examined in the limit of a large number of colors and flavors. The 
Hamiltonian matrix is given in a Fock space spanned by 't Hooft meson states and, for the case 
of zero fermion mass, a submatrix is diagonalized numerically to give the low-lying spectrum as a 
function of Nf/Nc- Pair creation effects generate bound states which are complicated mixtures of 
components of different meson number. There are a number of nontrivial zero modes; in the massive 
part of the spectrum some states tend to a well-defined Nf ^> Nc limit while others become unstable 
and disappear. The masses of most states remain remarkably constant over a large range of Nf/Nc- 



I. INTRODUCTION 



Since its inception, the large Nc expansion Q| || has afforded a host of valuable insights into QCD dynamics 
beyond standard perturbation theory in the coupling constant. Besides providing a solvable example of a nonabelian 
gauge theory in the case of one space dimension 0, it has been significant e.g. in clarifying the relation between the 
Skyrme model and QCD and in qualitatively explaining the Zweig rule ||. On the other hand, the number of colors 
Nc is not the only parameter of QCD which can be regarded as large. The same procedure may also be applied to 
, the number of flavors Np. In the real world, especially at higher energies, where the heavy quarks contribute to the 
^vq ' dynamics, the effective Np/Nc ratio becomes larger than unity. Thus it seems natural to investigate in what way the 
t-H , presence of a large number of flavors modifies the large Nc picture and how those parts of large Nc phenomenology 
which appear to be successes survive the modification. Early work in this direction was done by Veneziano and 
some more recent explorations into the field of multiflavor gauge theory have been carried out e.g. in connection with 
C| ■ bosonization treatments H || . 

fi i' Particularly the low-dimensional models such as QCDi + i with a large number of colors become considerably 
JL , simplified in the framework of light-cone quantization Q . At this point it should be remarked that the light-cone 
treatment corresponds in some respects to an effective theory with well-defined limitations, as was elucidated in Q. In 
the work presented here, no attention will be paid to such questions. The results will be obtained entirely within the 
usual naive light-cone framework and no firm claim is made as far as their validity in normal coordinates is concerned. 
The prevailing experience gives rise to the hope that the light-cone treatment does yield correct excitation spectra. 

The present foray into the multiflavor world focuses purely on the mesonic spectrum of the model with zero quark 
mass matrix. The Hamiltonian is calculated in a (truncated) basis of states made up of 't Hooft mesons, which 
diagonalizes the Hamiltonian for Np — 1, Nc — ► oo. Increasing the number of flavors such that Np/Nc remains finite 
as Nc — * oo induces nonvanishing interactions between the 't Hooft mesons; the evolution of the low-lying spectrum 
as a function of Np/Nc can be followed by rediagonalizing the Hamiltonian. Here, only the flavor singlet sector is 
discussed; however, the generalization to finite flavor density presents no conceptual problem. Since the numerical 
work up to now was carried out with limited computational resources, only the one-meson and two-meson (i.e. two- 
and four-parton) sectors could be included. 



II. THE MESONIC HAMILTONIAN 



The Hamiltonian of multiflavor QCD\+\ in the light-cone gauge with zero quark mass matrix reads in momentum 
space 
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Here, i and j are the color indices, a and b the flavor indices. Note that the fermion fields do not carry a Dirac 
index, since one of the components completely decouples in the light-cone gauge. The Hamiltonian has already been 
normal-ordered, by applying Wick's theorem, with respect to the vacuum of the theory, which is the perturbative one: 

9l»|0)=0forp<Q (2) 
q ai (p)\0) =0forp>0 (3) 

The singularities in the integrands finally are defined via the principal value prescription. As already mentioned in 
the introduction, all subleties associated with the singular nature of the light-cone description are disregarded here; 
the reader is referred to (]]]. 

The Hamiltonian (^) conserves baryon number, total momentum, and overall color. Here, the focus will be on 
the color singlet sector JTo| with zero baryon number and total momentum denoted by 2K. Due to the confining 
Coulomb potential, the physical spectrum consists entirely of color singlet quark- antiquark bound states, i.e. mesons, 
and conglomerates thereof (ll]] [Oj. It thus seems meaningful to consider the Hamiltonian matrix in a Fock space 
spanned by color singlet mesonic excitations above the (perturbative) vacuum. Mesons are created by the bilinear 
operators 

Ml iab = ^J q dpf n (Q, P ) iUp)m(p Q) (4) 

where the /„ are a complete set of bound state wave functions. A particularly suitable such set are 't Hooft's wave 
functions, 

f n (Q,p) = * ^(p/Q) (5) 

where the <p n satisfy the equation 
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with the boundary condition 



The mass scale too is given by 
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and [j.^ is the invariant mass squared of the n-th 't Hooft meson. Apart from the exact zero- mode solution <fio(x) = 1, 
a suitable orthonormal set in which to expand the 4> n in order to solve (((]) numerically is given in this case by cosine 
functions. 

With this choice, the Hamiltonian is already diagonal in the basis defined by the MQ nab in the case that Np 
remains finite and Nc — ► oo. Note that the operators M' are not suitable for a rigorous bosonization of the theory 
since they only obey canonical bosonic commutation relations up to additional terms of the order 1/Nc- It is now a 
straightforward, albeit lengthy, procedure to calculate the Hamiltonian matrix by commuting through quark operators. 
One obtains 

(2K, n, a, b\H\2K, m, c, d) = S ac S b dS nm ^ (9) 
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Here, in the kinetic part of the two- meson-two-meson matrix element, a convention of definite ordering of the momenta 
in the states was assumed, namely Q, Q' > 0. This avoids double-counting two-meson states. Also, in this kinetic 
part, l/iVc-suppressed terms in the normalization of the states were neglected. All interaction parts, on the other 
hand, are valid for any Q, Q' and no terms were neglected. 

The processes described by f mm ' n , fmm^nn'^. 3 ^^ fmm/nn' are associated with a flavor flow represented diagram- 
matically in graphs (A),(B) and (C) of figure (|l|), respectively. These graphs clarify the large Nf counting associated 
with each vertex; this will crucially influence which contributions survive in the limit Np/Nc constant, Nc — * oo. 

Note that the Hamiltonian can only change the number of mesons by one. To change it by two, it would have 
to contain a contribution consisting of four parton creation or annihilation operators with respect to the vacuum; 
this, however, would completely commute through any meson operators to the right or left, respectively ultimately 
annihilating the vacuum in relation to which the Hamiltonian is already normal-ordered. Equations (p|)-(|Tl|) completely 
define the Hamiltonian in the subspace containing up to two mesons. But for the fact that the meson operators M' 
only obey canonical commutation relations up to terms of order 1/Nq, the three- meson and four- meson interactions 
described by (||)-([ll]) would already represent the full dynamics in mesonic language. This is because the Coulomb 
interaction, as is clear from the original quark representation (Q), can only modify at most two out of an arbitrarily 
long string of meson operators as it is being commuted through to eventually annihilate the vacuum. However, the 
1 /A^c-suppressed additional terms in the meson commutators introduce additional interactions in larger clusters which 
in general must be treated on the same footing as the vertices derived above. Since the Hilbert space will here be 
truncated to the one- and two-meson sectors, these terms need not be considered further in the present calculation. 



III. THE FLAVOR-SINGLET SECTOR 



The above considerations still apply to an arbitrary structure in flavor space; the results will now be specialized to 
the flavor singlet sector. A flavor singlet meson can be constructed by the superposition 



\Q,n) 
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(13) 



whereas a flavor singlet two-meson state may consist either of two singlet mesons or of two mesons carrying flavor 
spin one coupled to a singlet: 



\Q, n; Q', n') ss = -j- ]T M^M^n' 

a,o 
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Carrying out the flavor algebra, the Hamiltonian matrix in these states takes the following form, where only terms 
which remain finite in the limit Np/Nc finite, Nq — > oo are kept: 

s (2K,n\H\2K,m) s = 6 nm ^ (16) 
s (2K,n\H\K+%,m;K-%,m')ss = (17) 
s{2K, n\H\K+%, m; K-%,m') TT = ^^^^(l + (-1)™+™'+"+*)/^ (18) 

(19) 
(20) 
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Here, the symmetry behavior of 't Hooft's wave functions <p n (x) around x — 1/2 was used, making explicit the 
selection rules for the meson excitation quantum numbers. Physically, this is the manifestation of CP invariance. 
These selection rules decouple states containing single-meson components of even excitation number from those 
containing single-meson components of odd excitation number. Note that this feature does not persist in the flavor 
non-singlet sector. 

The Hamiltonian matrix obtained above exhibits in very transparent form the physics of the large Np world: Flavor 
singlets do not interact among each other. Nonvanishing interactions only occur within flavor singlets, which resemble 
loops built out of alternating color strings and flavor "strings" in the sense that the two quarks at the end of a flavor 
string must be coupled to a singlet. Within these loops, color strings may break into two, creating a flavor string in 
between. Flavor strings are not necessarily binding; however, the diagonalization of the Hamiltonian will show that 
there are indeed bound states made up of more than one meson such that the flavor string effectively acts as a bond. 
Note the remarkable similarity of this picture to the model with one flavor, but adjoint color || ||. 

In the original quark Hamiltonian, flavor does not appear as a quantum number related to the interaction; thus, the 
notion of a flavor string introduced above does not carry physical meaning in the same sense as the color string, which 
in fact represents a region of nonvanishing chromoelectric field energy density. The flavor string "force" is generated 
entirely by the quantum fluctuations, which manifest themselves in the diverging polarizability of the vacuum as 
Np — ► oo. The weakness of the color forces between the mesons is offset by the ease with which quark- antiquark pairs 
are created. Formally, this is reflected in the fact that only graphs which imply the existence of free flavor indices to 
be summed over have a chance to compensate the usual l/iV<7-suppression of meson vertices by powers of Np in the 
numerator. These are exactly the graphs in which a (flavor singlet) quark-antiquark pair is created or annihilated, 
i.e. a color string is broken into two or vice versa. Thus, only the graphs (A) and (B) in figure (JjJ) survive; only the 
form factors / mm <„ and f^m'tin' en ter the Hamiltonian matrix (|l6|)-(|2l|). 

In contrast to the one-flavor theory, where pair creation is not enhanced, and which is simply solved by 't Hooft's 
mesons, the bound state spectrum will now consist of superpositions of states containing different numbers of mesons. 
In the case of the one-flavor, adjoint color theory || it was found that the low-lying spectrum is still very well 
described by a truncation to few quark-antiquark pairs. This gave rise to the hope that the restriction to one- and 
two-meson states imposed in the present calculation would also give an adequate description of the first few bound 
states. Ultimately, this point will have to be checked by enlarging the Hilbert space. 

Numerical evaluation and diagonalization of the Hamiltonian matrix (|l^)-(|2^) gives the low-lying spectrum as a 
function of Np/Nc- I n the actual calculation, the two-meson states were supplied with a relative momentum wave 
function, i.e. the physical basis used was of the form (cf. the convention introduced further above of definite ordering 
of momenta in the kets): 

r 2K ^ 

dQ ipi(Q/2K)—= QK+Q/2, to; K-Q/2, m') T T ± \K+Q/2, to'; K-Q/2, m) TT ) (22) 
V2 

where the ipi were chosen to be (properly normalized) cosines (with the plus sign) and sines (with the minus sign). 
Note that for m = to' only even (cosine) relative momentum wave functions are possible due to the Bose symmetry. 
Note also that the form factors ( |l2|) contain logarithmic divergences for zero quark mass which pose no problem if 
a continuous relative momentum wave function is introduced as above, whereas they make a direct discretization of 
the momentum difficult. The basis used in the computation consisted of eight one-meson and 26 two-meson states. 
Bound states were identified by subtracting the expectation value of the kinetic energy from the energy eigenvalue 
and checking for a negative result. 
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IV. NUMERICAL RESULTS 



The numerical calculation yielded four nontrivial zero modes besides the zeroth 't Hooft meson \2K, 0)<j, which is 
an eigenstate of the Hamiltonian for all Np/Nc, as can be easily seen from the three- meson form factor f m m'n (cf. 
(|l2])). Furthermore, a spectrum of massive bound states, plotted in figure (||) as a function of Np/Nc, was found. 
Some of the massive states disappear for intermediate values of Np/Nc and reappear with very similar energy and 
mesonic composition. In such cases the range where the state disappears is bridged in figure (J2J) by the dotted lines 
to indicate the kinship between the spectral trajectories in different regions. 

The different classes of states obtained can each be described in more detail as follows. The nontrivial zero modes 
are of the structure 

1 r 2K 

V a„|2A» s + ^= / dQ X (Q/2K)\K + Q/2, 0; K - Q/2, 0)tt (23) 
nodd ^ K J o 

where the restriction to odd n comes simply from the selection rules already discussed above in connection with 
equations (|n])-(|2l|). It should be noted that the Hilbert space used in the present calculation contained the first four 
relative momentum states made up of two ground state 't Hooft mesons and also the first four single-meson states of 
odd excitation number. In this basis, four nontrivial zero modes were found. It thus seems natural to expect that the 
full theory in fact contains an infinity of zero modes of the form (p3|). For small Np/Nc, the one-particle admixture 
vanishes as Np/Nc- Such a scaling behavior is to be expected; for Np/Nc ~ > 0, the Hamiltonian becomes diagonal 
in the mesonic basis and thus the eigenstates are dominated by a single (either one- or two-particle) Fock state, with 
./Vp/iVc-suppressed admixtures obtainable in perturbation theory. The same behavior is exhibited by the binding 
energies. This feature of course also arises in the massive part of the spectrum to be discussed further below. 

At Np = Nc, the one-particle content ranges from 40% for states in which these one-particle admixtures are made 
up of low-lying 't Hooft mesons, downwards as more excited mesons participate. The binding energies at Np/Nc are 
exactly 2m^/K (up to the numerical accuracy of the calculation) for three of the zero modes, rising to 2.hm} ) /K for the 
fourth. In the opposite limit of large Np/Nc on the other hand, the two-meson admixture is suppressed as Nc/Np 
and the states become mixtures of odd single 't Hooft meson states. As before, this scaling is easy to understand: 
For large Nf /Nc, the two-meson-two- meson interaction Hamiltonian, which is positive (this will be discussed further 
below), becomes dominant, and thus it is not favorable for a state to contain a high two-meson proportion. In order 
to maintain a meaningful balance between the different parts of the Hamiltonian, the two-meson amplitude must 
be proportional to y/ Nc /Np , i.e. the two-meson content of the states must scale as Nc/Np. Again, this property 
of course persists for the massive part of the spectrum. The binding energies of the zero modes must accordingly 
approach constants. Numerically one obtains two states with a binding energy of Sttiq/K for large Np/Nc, one with 
binding energy 2Qrn^/ K, and one with 2§m^ j /K. It has not been possible to give analytical solutions for the nontrivial 
zero modes; there must be subtle dynamical mechanism responsible for suppressing transitions from the single-meson 
admixtures in (^) to the massive part of the spectrum. It should be noted that some of the massive bound states to 
be discussed below also contain sizeable admixtures of odd single-meson states. 

Continuing on to the massive part of the spectrum, the states containing odd single-meson components appear as 
the third and the fourth massive state in the Np/Nc "C 1 region (cf. figure (^)). In detail, these two states however 
exhibit some differences: The fourth state only exists at low Nf/Nc, and in the region shown consists of about equal 
parts of 1 5) and |1, 1) (the quantum numbers in the kets merely indicate the 't Hooft meson content, whereas the 
momentum quantum numbers have been suppressed). For very low Np/Nc, though, the state becomes dominated by 
the |1, 1) part. By contrast, the third massive state becomes dominantly a one-meson state in the Np/Nc ~ * limit, 
namely a |3). As Np/Nc rises, it acquires admixtures of |0, 0), |1, 1) and |0, 2) so that the |3)-component is still found 
with 55% probability at Np = Nc- Whereas all other states found are dominated by two to three Fock components 
even in the transition region around Np — Nc, this state becomes a strong mixture. This is especially pronounced 
where the state becomes unstable. Remarkably, as Nf/Nc rises beyond unity, the two-meson content rises further to 
reach 96% at Nf/Nc = 50; one third of this is the |0, 0) admixture. In accordance with the argument given further 
above, this ultimately upsets the balance between potential and kinetic energies. Thus, the state is very weakly bound 
above Np/Nc = 10 and finally disappears around Np/Nc = 75 instead of increasing its one-particle content to give 
a well-defined Np/Nc ~^ °o limit. The large range over which this state survives is remarkable, though. 

On the other hand, the rest of the massive states shown in figure (||) contain even single-meson states. Here again, 
one finds states which become dominated by a two-meson Fock component for Np/Nc - * (namely the second 
massive state, which becomes predominantly a |0, 1)) and others dominated by a one-meson component in this limit: 
The first massive state becomes a |2), the fifth an |8). The one-particle fraction of the first massive state falls to 
about two thirds at Nf = Nc due to admixtures of |0, 1) and then rises again to give a sensible Nf/Nc — > oo limit. 
Remarkably, in this region, its mass falls to zero. Apart from this last feature, the fifth massive state behaves very 
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similarly, acquiring a |1,2) admixture of about one third in an intermediate regime of Np/Nc and then becoming 
again a pure |8). The second massive state acquires a steadily rising admixture of |4) (about 25% at Np = Nq) to 
rapidly become pure |4) for large Np/Nc- The state which is predominantly a |6) for Np/Nc — > 00 only appears 
at Np/Nc > 10; it is the fourth massive state in this region (cf. figure (|J)). This behavior is most probably due to 
the truncation of Hilbert space, i.e. the state may be dominated at lower Np/Nc by a two-meson component not 
included in the calculation. 

Summing up, the gross features of the spectrum can be characterized as follows: On the one hand, there are the 
zero modes, which consist predominantly of |0,0) with growing admixtures of odd single-meson states as Np/Nc 
rises until these admixtures become dominant for large Np/Nc- Concomitantly, the massive states containing odd 
single-meson components become unstable as Np/Nc is increased, which is signalled by these states becoming strong 
mixtures without a dominant Fock component. On the other hand, the massive states containing even single-meson 
admixtures display a well-defined large Np/Nc limit, where they are dominated by just these even single-meson 
states. The lowest of these states becomes massless for Np/Nc ~ > 00 ■ The reason for this exceptional behavior is 
unclear; all other states display a remarkably stable mass over a wide range of Np/Nc- Also this stability of the 
masses is a puzzle; the binding energies of the states on their own vary quite strongly, as is evident from the values 
indicated in figure (0). 



V. COMMENTS AND OUTLOOK 



Some comments are in order regarding the results presented above. The first of these concerns the low Np/Nc 
limit. The reader will have noticed that not all of 't Hooft's mesons appear in the bound state spectrum in this limit. 
To be precise, only |0), |2), |3) and |8) are present; |1), |4), |5), |6) and |7) are missing. These states of course do 
appear in the numerical diagonalization, with small admixtures of two-meson states, but with a positive expectation 
value of the potential energy. Thus they are not bound and must rather be viewed as resonances. For very small 
Nf/Nc, this distinction becomes meaningless. An experimenter in 1+1 dimensions, working in a finite laboratory 
for a finite time, will observe all of 't Hooft's mesons as Np/Nc ~> 0. Strictly speaking, however, if for fixed Np/Nc, 
the experimenter waits long enough, he/she will only find some of these mesons as stable particles. In the case of 
the higher mesons, the picture may be strongly influenced by the truncation of Hilbert space, but the absence of the 
|1) at least seems conspicuous. On the other hand, the bound states become very weakly bound for Np/Nc ~ ¥ 0j so 
that the experimenter would indeed observe the |0, 1) (this is the second massive state) as two essentially independent 
ground state and first excited 't Hooft mesons, respectively, in a finite experiment. 

As a further comment, it should be stressed that, to generate bound states in the model, it is crucial to combine 
one-meson and two-meson states. Diagonalizing the two-meson sector on its own does not give binding, i.e. the 
two-meson-two-meson interaction matrix is positive, as mentioned above. This seems plausible, for the following 
reason: If there was a pure two-meson state with a negative expectation value of the interaction, then one would have 
a state with arbitrarily negative overall energy by simply choosing Np/Nc suitably large such that the interaction 
dominates the kinetic part. On the other hand, the ground state (vacuum) energy of the theory has been defined 
to be zero, so that a state with negative energy should not be possible. Conversely, including one-meson states can 
very well generate binding through the three-meson vertex, and the energy is stabilized to remain positive by the 
two- meson-two-meson vertex; this is indeed precisely what happens, as mentioned above in the detailed discussion of 
the spectrum at large Np/Nc- If one writes a mesonic Hamiltonian in terms of the operators M' and their adjoints 
such that it correctly reproduces the Hamiltonian matrix (|9|)-([l"l|), it is the cubic term that may become negative, 
and the quartic term which stabilizes the theory. 

Summing up, the large Nc, large Np theory contains a considerably wealthier phenomenology than the one-flavor 
theory. It represents a more sophisticated model for mesons, incorporating pair creation effects and thus mimicking 
the presence of a meson cloud, thought to be an important ingredient of three-dimensional hadrons. Despite the 
complex behavior of the model, many of its features seem to be susceptible to numerical analysis, given the power of 
modern-day computers. In this respect, a significant improvement on the treatment presented here should be possible, 
since the latter only made use of two workstations. By using a mesonic basis, which already incorporates a sizeable 
part of the dynamics, the bulk of the computation lies in the evaluation of the Hamiltonian matrix rather than its 
diagonalization. The matrix elements, however, can be evaluated completely independent of each other, making a 
large improvement in speed possible by simply using parallel processors. Including sectors with a higher number of 
mesons demands the computation of higher-dimensional integrals; this will only be feasible beyond a certain point by 
employing Monte Carlo methods. 

Besides improvements of the numerical aspects, there are a number of other extensions of this work which seem to 
merit attention. It would be interesting to understand in more detail the dynamics leading to binding. This endeavor 



G 



would profit greatly from some analytical understanding of the nontrivial zero modes. Especially the zero mode 
sector may to some extent be tractable analytically, at least in some approximation, since the form factors simplify 
considerably if some of the participating mesons are ground-state ones. 

Furthermore, the approach followed here becomes the most questionable in the high Np/Nc limit. There, the 
coupling between sectors of different particle number is the strongest and the inclusion of Fock components with three 
or more mesons could substantially alter the results. It would be helpful to devise an approximation scheme which is 
particularly suited for this potentially interesting new limit, i.e. one which makes explicit use of the smallness of the 
parameter Nc/Np. 

With a view to approaching more realistic theories, it may be helpful to work out the perturbative corrections in 
1/Nc- Finally, note also that the model is easily generalized to non-singlet flavor (flavor spin one objects will resemble 
open strings) as well as non-zero quark masses. In fact, in the latter case, the Hamiltonian matrix retains exactly 
the form (|9|)-(|lT|); only the 't Hooft wave functions entering the form factors and the corresponding meson masses are 
modified. It is to be expected that quark masses suppress pair creation related effects. Thus, the physical spectrum 
should remain closer to the 't Hooft one for given Np/Nc- The zero quark mass model considered in this work is the 
case furthest removed from the well-studied one-flavor theory, since it maximally allows for the pair creation processes 
which represent the essence of the model's nontrivial dynamics. 
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FIG. 1. Flavor flow associated with the form factors / mm ' n , /mm'nn' an d fmrrJ~nn'> respectively. 

FIG. 2. Invariant square masses of the massive bound states as a function of Nf/Nc- Dotted lines are merely to guide the 
eye in connecting spectral trajectories of states which disappear at intermediate values of Np/Nc- The figures indicate the 
binding energies of the states at selected points in the same units, i.e. K{H po t) /2mJ. 
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